Gliding dislocations in a driven vortex lattice 
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The dynamics of dislocations in a two-dimensional vortex lattice is studied in the presence of a 
pinning potential and a transport current. In a vortex lattice drifting with velocity v a glide velocity 
V<j of the dislocation with respect to the vortex lattice is found to decay like Va ~ v~ 4 for large 
drive. From this result the velocity for the crossover between a regime of coherent elastic motion 
and a regime of incoherent plastic motion of vortices is estimated. 

PACS numbers: 74.60.Ge, 61.72.Bb 
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Periodic structures subject to quenched disorder are a 
generic model for a vast number of physical phenomena. 
The recent upraise of interest was motivated by the dis- 
covery of high temperature superconductors (HTS) that 
set a physical realization of the random system with tun- 
able parameters and which are easily accessible in ex- 
periments |lj] . As the statistical mechanics of the static 
periodic media in the presence of disorder reaches its ma- 
ture phase, the physics of non-equilibrium processes pro- 
poses a diversity of unresolved challenging questions. In 
the static case it is established that quenched pinning 
destroys the crystalline order of the vortex lattice on in- 
termediate spatial scales jl|-[|. The long-range order of 
the vortex lattice is reduced to the quasi-long-range order 
of the "Bragg glass" Mm. It is important to note that 
in the static case and for weak pinning, the topological 
structure of the lattice is preserved in a certain range of 
the phase diagram ||]. 

One of the central issues of non- equilibrium statisti- 
cal physics of dirty periodic media is the very nature 
of driven steady states. It was observed |j]-^| that the 
depinning transition can be accompanied by plastic de- 
formations of the driven periodic structures leading to 
the destruction of the moving solid. For driven vor- 
tex lattices it was proposed || that at larger drives a 
true dynamic phase transition occurs from an incoher- 
ently moving non-equilibrium state to a coherently mov- 
ing solid phase. This conclusion was later extended to 
driven charge density waves [ |10| . Such dynamic transi- 
tions have indeed been observed in experiments pd| and 
simulations |],|l2|]. Note, that even in the limit of strong 
driving disorder severely affects the lattice structure on 
large length scales (in dimensions d < 3), which becomes 
strongly anisotropic and keeps glassy features |13|H15|. 

A phenomenology of this dynamic melting is built on 
the notion that the driven periodic medium experiences 
pinning as temporally fluctuating distortion in the co- 
moving frame. A convenient heuristic tool to investi- 
gate the dynamic transition was introduced in Ref. ||, 
where the effect of temporarily fluctuating disorder was 
described as an additional thermal noise with an effec- 
tive "shaking temperature" . This temperature increases 



when the current is reduced and can lead to a "dynamic 
melting" of the lattice structure. While the phenomeno- 
logical approach enables to estimate the position of the 
transition line on the phase diagram, the nature of the 
dynamic melting remains an open question. A micro- 
scopic understanding of dynamic melting requires the 
analysis of topological defects, which become free at the 
transition. In this Letter we present a first step in this 
direction focusing on the study the dynamics of inde- 
pendent dislocations in the driven two-dimensional vor- 
tex lattice. Note that in a homogeneously driven vortex 
lattice without pinning dislocations and vortices move 
with the same average velocity. Inhomogeneous pinning 
of different parts of the vortex lattice causes shear strains 
which generate a fluctuating force acting on dislocations. 
A motion of dislocation relative to the vortices arises from 
correlations of this force on intermediate scales. We ex- 
amine the glide motion of the dislocations in the driven 
vortex lattice in the limit of large vortex drive, weak dis- 
order, and well below the melting temperature of the pure 
vortex lattice. We calculate the dislocation velocity per- 
turbatively and find that dislocations are retarded with 
respect to the vortices. The velocity where dislocations 
get pinned by disorder provides a non-equilibrium esti- 
mate for the crossover between elastic and plastic vortex 
flow. 

In order to derive the equation of motion for the dislo- 
cation we first consider the dynamic response of the vor- 
tex lattice to disorder and thermal fluctuations. We re- 
strict our discussion for simplicity to the two-dimensional 
vortex lattice. This lattice is supposed to move with an 
average velocity v as a consequence of the Lorentz force 
i"l generated by electric transport currents. We label vor- 
tices by their perfect lattice position R in the comoving 
frame. Their actual position in the laboratory frame is 
denoted by r. Then a displacement can be defined as 
usual by u(R, t) = r(R, t) — vt. 

Neglecting nonlocal effects the elastic energy of the 
vortex lattice reads 

Hc ^\JWf {^l^l 2 + c ^ 2 |u T (q)| 2 } , (1) 
where q runs only over the first Brillouin zone and u L,T 



are the longitudinal / transverse projections of the dis- 
placement 

field u(q,w) = J dt Oo/ B ) Er e™*-'*™ u(R,t) and 
4>o/B is the area per vortex. The pinning potential V(r) 
is supposed to have correlations V(r)V(r') — A(r — r'). 
The pinning energy of the lattice is H P i n — J2n ^( r (R))- 
The motion of vortices follows the overdamped 
Langevin equation 



7?r(R,t) = - 



d 



dr(R,t) 
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H pin ] + f L + £(R, t) (2) 



with the Bardeen-Stephen friction coefficient 77 and 
a thermal noise correlation (£ a (R, t)£p(R', t')) = 
2r)T6 a0 6 RW 5(t-t'). 

In the limit of large driving forces fn, the response of the 
vortices to the pinning potential vanishes. Then one has 
to leading order v = ih/v and the external force acting 
on the vortices is f(R,t) := £(R,i) - W(R + v£). Its 
correlations read 



</ Q (q, «)/^(q', «')> = Sa/J(q, w)*(q + q')*(w + w'), (3) 



S Q/ j(q, w) = 2 V T^S aP + J2 k a k p A(k)S(uj + k • v) 

Q 



where k := Q + q with a reciprocal lattice vector Q. In 
linear response H the vortex displacements are obtained 
from the forces f through the response function with 
longitudinal part r L (q, u>) = [— irju + cn(<j)o/B)q 2 ]~ 1 



and transverse part T (q, uj) = [— irjui + CQ6(tfio/B)q 2 ]~ 1 . 
Typically the longitudinal response is small compared to 
the transverse response since en > cqq. 

As one can immediately recognize from expression ( H) , 
pinning forces are more relevant than thermal forces fl3f | 
on large length and time scales: thermal fluctuations give 
a finite contribution to 2 Q ^(q = 0, uj = 0), whereas pin- 
ning forces give singular contributions for Q ^ with 
Q • v = 0. This effect is due to the periodic structure of 
the vortex density and leads to the destruction of transi- 
tional long-range order for all velocities v in dimensions 
d < 3 ]13||L5|]. For any finite driving the correct asymp- 
totic behavior of the lattice roughness cannot be obtained 
in a linear response approach. However, as we show be- 
low, dislocation dynamics is governed by the elastic re- 
sponse of the vortex lattice on intermediate scales where 
this approach is sufficient. 

The topological structure of the vortex lattice is ex- 
pected to be destroyed by dislocations and disclinations, 
at least for sufficiently high temperature, strong disorder, 
and/or weak driving. Within an elastic approach vortices 
are found to move like in channels (lqJl^l when the lat- 
tice is moving along one of its principal directions. The 
displacement field component having the strongest fluc- 
tuations is the component parallel to the velocity pT].|n>l ■ 
This favours dislocations with Burgers vectors parallel to 
the velocity. 

The motion of such a dislocation is easiest in the glide- 
direction parallel to the Burgers vector b with length of 



the lattice spacing a w \JWqJB. In this direction it can 
move by slips of individual vortices. A motion of the 
dislocation perpendicular to the Burgers vector (climb) 
is possible only in combination with the creation of va- 
cancies or interstitials. In other terms, glide is controlled 
by the Peierls barrier which is much lower than the en- 
ergy barriers for dislocation creep. In what follows we 
thus can ignore dislocation climb. The dislocation mo- 
tion is then one-dimensional in the direction parallel to 
the Burgers vector and to the velocity. Its position along 
this direction in the moving vortex frame is denoted by 
X(t), and we can consider Y(i) = fixed. 

According to the theory of dislocations MM internal 
stresses 



cr a p = 2c 66 y a ui3 + (en - 2c 66 )V 7 u 7 <5 Q 



(4) 



of the lattice, due to thermal fluctuations and pinning, 
induce the so-called Peach-Koehler (PK) force 



K a — Sapo-p^b^ 



(5) 



acting on the dislocation, where e denotes the antisym- 
metric tensor. In addition, the dislocation is subject to 
a periodic potential W(X) w — (W p /2) cos(27rX/a) with 
amplitude W p w \Q~ 2 cqqo? representing the Peierls bar- 
rier for glide motion. The resulting equation of motion 
for the dislocation reads 



77 d X(i) = K x (X(t),t) - V x W(X(t)). 



(6) 



Here the dislocation friction coefficient is taken equal 
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to its low velocity value rjd rs ( 77/871- ) In (l/nda 2 
for an average dislocation density n^. We assume that 
the vortex lattice is still ordered on intermediate length 
scales, i.e. the density of dislocations is small compared 
to the vortex density and ?7d 3> 77. 

The fluctuating PK force includes thermal and disorder 
induced parts. The corresponding spatio-temporal cor- 
relator follows within linear response straightforwardly 
from Eq. (g) and has has a rather complicated form. For 
the further analysis we single out the most important 
contributions. 

One contribution arises from the thermal noise £ act- 
ing on the vortices and it is present even in the ab- 
sence of disorder and a driving force. It can be de- 
scribed as effective thermal noise £d(£) with correlations 
(£d(t)£d(t')) ~ 2r]dTS(t — t') experienced by dislocations. 

The other contribution comes from pinning forces gen- 
erating fluctuating displacements u(R, t). As pointed out 
in Ref. [O] , the modes Q with Q • v = of the force cor- 
relator (H) generate displacements, which are stationary 
in the laboratory frame, u(°)(R,i) = u.(°)(R + vt). The 
associated contribution to the PK force is also stationary 
in the laboratory frame and can be the only origin of a 
pinning of dislocations by disorder. The correlations of 
these contributions K x at Y = are 



Ki°\q x )Ki°\q' x ) 



x(<lx)S{q x + q' x 

16a 2 c§ 6 A (r)V\(/,. 
iitarfv 2 



(7) 
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Here A := YIqq =o Qy^iQ) an d ^he obtained formula 
is valid for i]v 5g> c\\a. Although the mode Q = also 
gives a stationary contribution, it is negligible as com- 
pared to the modes Q ^ with Q • v — 0, because of 
en S> C66- We ignored the modes Q • v ^ which repre- 
sent a shaking of the dislocation due to the discreteness 
of the vortex lattice and are expected only to renormal- 
ize the effective dislocation temperature. This effect is 
negligible for large v, and does not contribute to the dis- 
location glide. 

Having made the above approximations we rewrite the 
equation of motion (g) determining the dislocation glide 
velocity in a form 



Vd 



X = K^(X + vt) 



ttW„ 



2nX 



■& 



(8) 



We restrict ourselves to temperatures T <C W p , i.e. well 
below the melting temperature Tkt ~ (3c66a 2 )/87r of 
the pure lattice J20|. In this case the dislocation moves 
via thermal activation from one minimum of the Peierls 
potential to a neighboring one favored by the PK force. 
In the driven lattice disorder biases this motion in the 
direction opposite to the vortex drift. 

For weak disorder A <C T 2 this effect is conveniently 
calculated perturbatively in the dynamical functional for- 
malism pl| , where the dislocation position X and its re- 
sponse field X are assigned a statistical weight exp(— .4) 
with a dynamic action 



A = Ao + A u 



(9) 



A 



dt{ m TX 2 + iX[r) d X + VxW(X)]}, 
,'!', I dt'XX' X {X - X' + v(t - t')). 



Here the disorder average is already performed and ab- 
breviations X = X(t), X' = X(t') etc. are used. 

In the absence of disorder, dynamics is controlled by 
the free action Aq. This action is minimized by instanton 
trajectories which start from a minimum/maximum of W 
at t = — oo and end at a neighbored maximum/minimum 
at t = oo: 



X\ (t) = a 2 — arctan 



cxj) | a 1 - 



Xi(t) = - 



2r) d T 



[ThXi(f)+^xW(Xi)], 



(10a) 
(10b) 



where t\ — a 2 r]d/2iT 2 W-p is an instanton time scale and 
cti,(T2 = ±1 specifies the type. For o\ = ±1 motion is 
up or down from a minimum to a maximum of W or vice 
versa. For <j\&i = ±1 motion is forward/backward in X 



direction. Both "down'-instantons have X^ own (t) 
and thus „4$ own = 0. Both "up"-instantons have 



^ p = W p /T, 



X?»(t) = 



2-rrtiT cosh(i/ii) ' 





(11a) 
(lib) 



Eq. (JT^a) gives the usual statistical weight of thermally 
activated instantons. 

Now one can easily calculate the disorder corrections 
to the instanton actions in perturbation theory. To eval- 
uate A\ in the lowest order in A one has to use the un- 
perturbed instanton solutions and recognizes immediately 
that Af own = since A >down = 0. Near its saddle-point 
the dislocation has a velocity V[ = v + Xi(0) relative to 
the disorder, which is different for "forward" and "back- 
ward" instantons. In terms of this velocity the action 
correction reads 

Al ™\ I dt fdt'XX'x(vi(t-t')) 



2A [2a 2 c m 



3n 2 T 2 yirnvviti 



3/2 



(12) 



The resulting dislocation glide velocity Vd '■= 
(X) is obtained from the difference of activation 
rates for forward and backward instantons, V^ ~ 
(a/£ )[exp(-.4 fw ' u P) - exp(-.4 bw ' u P)]. The attempt fre- 
quency of the dislocation t rs a 2 ru/irW p is estimated 
according to Kramers rate theory [ B2[ . In terms of the 
dislocation diffusion constant Dq = (a 2 /to) exp(— *4o) m 
the pure lattice we eventually find a drift velocity 



Vd 



L> A / 2a 2 c 66 



vt\T 2 \ir 3 r]v ti 



3/2 



(13) 



in the limit v — > oo, where Ai <C 1. 

Several properties of this resulting drift velocity are 
noteworthy: it decays with the fourth power of veloc- 
ity, which is much faster than Vd ~ v^ 1 that one might 
expect on the basis of symmetry arguments. The drift 
velocity is independent of Cn for tu » cqq, where only 
shear modes of the lattice get excited. In the consid- 
ered temperature range T <S^ W p the drift velocity arises 
from a small bias of instantonic hops. Since this pro- 
cess is thermally activated the drift velocity increases 
monotonously with temperature. The asymptotic form 
(O) holds only down to a certain minimum temperature 
determined by the condition A\ *C 1, which however is 
satisfied for large velocity. 

From Eq. (H3) we can estimate the critical velocity 
where a crossover from "elastic" to "plastic" flow takes 
place. In the "elastic" flow regime vortices flow essen- 
tially coherently. Topological defects may be present due 
to thermal fluctuations. Although we found that dislo- 
cations are always retarded with respect to the vortices, 
this effect is small for high drift velocities v 3> |Vd|. On 
the other hand there is the "plastic" flow regime, where 



dislocations are essentially pinned by impurities. Such 
dislocations are a source of rearrangements of neighbor- 
hoods between vortices. This plasticity becomes signif- 
icant for 1 — |Vd|/w <C 1. Although our derivation of 
Eq. (|l^) is strictly valid only for large velocities, we ex- 
pect it to give a qualitatively correct estimate for the 
velocity 



D q Aq (2a 2 c 66 



t{T 2 \ TrVi 



3/2' 



1/5 



(14) 



of the crossover between both regimes. 

Although the instantaneous forces acting on the dislo- 
cation arise from local stresses in the vortex lattice, the 
dislocation glide velocity results from the correlations of 
these stresses on intermediate scales \X\ < vt\ for large 
drift velocities: According to Eq. (|ll|b) the quantity 
Xi{t) decays exponentially on the instanton time scale t\. 
Thus the main contributions to A\ in Eq. (||) come from 
x(A) for |X| < vt\. This scale increases ex v 1 , whereas 
the dynamic Larkin length, up to where the perturbative 
calculation of x is reliable, increases oc v 3 13 la]. It is 



for this reason that the calculation of x hr linear response 
approximation for the vortex lattice is sufficient for the 
dislocation dynamics at large drive. 

Relating our result v c oc A with the finding of 
a dynamic freezing velocity Vf oc Ao, one has to distin- 
guish the case of weak and strong pinning. For strong 
pinning v c < Vf. Expression (UJ) for v c is then formally 
invalid since the lattice is no longer ordered at these ve- 
locities. In this case we can describe only the regime 
v > Vf, where the vortex lattice preserves its crystalline 
order on large scales and dislocations move with the vor- 
tices. In the limit of weak pinning both velocities Vf < v c 
become small and the perturbative findings are no longer 
quantitatively reliable. Nevertheless, the smallness of the 
velocities indicates that the lattice is well ordered and dis- 
locations move essentially with the lattice at practically 
all drift velocities. Qualitatively, our finding suggests a 
crossover between practically elastic flow for v > v c and 
strongly plastic flow for v < v c . This crossover can be 
viewed as precursor of the pinning transition. 

To sum up, our result Vd oc t>~ 4 implies that the motion 
of the vortex lattice in disorder is always accompanied by 
plasticity, even if it is extremely weak at high drives. For 
weak pinning we have estimated the characteristic veloc- 
ity of the crossover between "elastic" and "plastic" flow 
of the vortex lattice. We expect these results to be valid 
also after the inclusion of interactions between disloca- 
tions. For example a dislocation-antidislocation pair will 
be retarded with respect to the vortex lattice essentially 
as a whole, since the resulting glide velocity is indepen- 
dent of the dislocation type. But even if this pair remains 
bound it can create point defects during its motion. In- 
teractions between dislocations are clearly relevant for a 
study of the qualitative nature of the dynamic melting 



transition Q , we will address this issue in a forthcoming 
publication. 
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